The relation between TMDs and PDFs in the covariant parton model approach 
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We derive relations between transverse momentum dependent distribution functions (TMDs) and 
the usual parton distribution functions (PDFs) in the 3D covariant parton model, which follow from 
Lorentz invariance and the assumption of a rotationally symmetric distribution of parton momenta 
in the nucleon rest frame. Using the known PDFs fl{x) and g\{x) as input we predict the x- and 
PT-dependence of all twist-2 T-even TMDs. 
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I. INTRODUCTION 

o 

^ 1-1 1-1 

TMDs [TJ [2] open a new way to a more complete understanding of the quark-gluon structure of the nucleon. 
fSJ Indeed, some experimental observations can hardly be explained without a more accurate and realistic 3D picture of 
Q the nucleon, which naturally includes transverse motion. The azimuthal asymmetry in the distribution of hadrons 
produced in deep-inelastic lepton- nucleon scattering (DIS), known as the Cahn effect ^, is a classical example. The 
intrinsic (transversal) parton motion is also crucial for the explanation of some spin effects jlHIB]. 

In previous studies we discussed the covariant parton model, which is based on the 3D picture of parton momenta 
with rotational symmetry in the nucleon rest frame [17.-. 2 6 . An important feature of this approach is the implication 
of relations among various PDFs, such as the Wandzura-Wilczek approximation between g\{x) and ^2(^) which was 
proven in [IF together with some other sum rules. Assuming SU{^) symmetry (in addition to Lorentz invariance 
and rotational symmetry) relations between the polarized and unpolarized structure functions were found [19 , which 
agree very well with data. In [20] transversity was studied in the framework of this model and a relation between 
Q^ transversity h\{x) and helicity gi{x) was obtained. In a next step we generalized the model to the description of 
OJ TMDs. We derived a relation between the pretzelosity distribution /i^^, transversity and helicity [22]. Finally, with 
the same model we studied all time-reversal even (T-even) TMDs and derived a set of relations among them [23] . 
Moreover, it was also shown that the 3D picture of parton momenta inside the nucleon provides a basis for a consistent 
^-H description of quark orbital angular momentum [21] , which is related to pretzelosity [26] . It should be remarked that 
^ some of the relations among different TMDs were found (sometimes before) also in other models [27-34 . 

The comparison of the obtained relations and predictions with experimental data is very important and interesting 
from phenomenological point of view. It allows us to judge to which extent the experimental observation can be 
interpreted in terms of simplified, intuitive notions. The obtained picture of the nucleon can be a useful supplement 
. to the exact but more complicated theory of the nucleon structure based on QCD. For example, the covariant parton 
C^l model can be a useful tool for separating effects of QCD from effects of relativistic kinematics. 

In this paper we further develop and broadly extend our studies ^24n25j of the relations between TMDs and PDFs. 
The formulation of the model in terms of the light-cone formalism [23] allows us to compute the leading-twist TMDs 
• • by means of the light-front correlators (/)(x, PT)ij [2^ as: 



(1) 



2 tr [7+75(?!>(a;,pT)] = S'Lfi'^la;, pr) + — ^5i^T K Pt), (2) 
I tr [ia^^j.Hx, Pt)] = hUx, pt) + S,^ hi£{x, pr) + ^^tP't - 1^¥'')S't ^^a^^^ ^ p^^^^^^ 

The main goal of this work is to derive relations between these TMDs and the usual PDFs fi{x) and gi{x). Similar 
tasks were recently addressed also in other approaches [35[[36^. The Section lllj is devoted to the case of the unpolarized 



TMD, and in Section III we discuss polarized TMDs. We use the obtained relations to calculate and discuss the 



numerical predictions for TMDs. The Section IV is devoted to our concluding remarks. 
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II. THE UNPOLARIZED TMD 

The distribution /i(x,Pt) is given in the covariant parton model approach by the expression 



f!ix,pT) = xM 



fdp^ 
J 



M 



x] =M G{f) . 



In the final step of (jZJ) we performed the p-*- -integration by rewriting the (5-function as 



(4) 



(5) 



The remarkable feature of the present parton model approach is that one can predict unambiguously the x- and 
PT-dependence of TMDs from the x-dependence of the corresponding ("integrated") parton distribution functions. 
The deeper reason for that is the equal (3D-symmetric in the nucleon rest frame) description of longitudinal (i.e. 

-dependence) and transverse (pT-dependence) parton momenta. Let us remark, that the invariant parameter x 
(Bjorken x) is tightly connected to both longitudinal and transverse parton through the ^-function. In the nucleon 
rest frame transverse momenta play for x an important role according to ([5|. 

The 3D momentum distribution G^{p^) was expressed in terms of fi{x) in previous works [18l[2T]. In order to make 
this work self-contained we present here an independent derivation. We start from the model expression for fi{x) 
which follows from the first equality in Eq. Q. For the remainder of this section we set the parton mass m ^ 0. 
Besides being a reasonable approximation, this step gre atly simp lifies the calculation though finite m-effect can be 
included [19 . Notice that if m is neglected then = ^Jp^ + Pt = ^- ^ow, instead of integrating over p^ as we did 
m Eq. (H, it is convenient to use spherical coordinates, and define the angles such that p^ = p cos 0, i.e. 



r^a/ \ c f P ~ P COS 



xMj'-fG^ip)8^ M 

/>27r nOO pi 

xM d(t) dppG\p) 

Jo Jo J-1 
r 

271 xM'^ 

Jo 



d cos 5 



p — p cos f 
M 



dpG^(p) ^[p--xM 



(6) 



The 0-function emerges because the integral over the (5- function obviously yields a non-zero result only if | cos^| < 1, 
and implies a lower limit for p-integral. Notice that there is also an upper limit, namely p < ^M, related to the 
fact that X < 1 [18 . This upper limit is natural in the covariant parton model in the nucleon rest frame because 
an on-shell parton can carry at most the momentum Pmax = \^ which must be compensated by all other partons 
going in the opposite direction, such that the center of mass of the nucleon remains at rest. Since this is an unlikely 
constellation the momentum distribution G^{p) vanishes as p ^ Pmax similarly as fi{x) drops to zero with x ^ 1. 
Thus we obtain 



: M 



= 2tiM^ j dpG^(p) 

h xM 



(7) 



and we reproduce the identity [TSl |2T 



d 

dx 





= -ttM^ G" ( 


X 





xM 



This result inserted in Q enables us to predict uniquely ffix^pr) from fi{x) as follows 



/i(^,Pt) 



1 



7rM2 dy 

with the dependence on x, pt given through the variable 



f!iy) 



(8) 
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FIG. 1: The TMDs /i(x,pt) for u- {upper panel) and d-quarks {lower panel). Left panel: /i(x,pt) as function of x for 
Pt/M — 0.10 (dashed), 0.13 (dotted), 0.20 (dashed-dotted hne). The sohd hne corresponds to the input distribution fl{x). 
Right panel: /i^(x,pt) as function oipr/M for x = 0.15 (soUd), 0.18 (dashed), 0.22 (dotted), 0.30 (dashed-dot Une). 
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FIG. 2: /i^(x,pt) as function of {pr/M)'^ for x = 0.15 (sohd), 0.18 (dashed), 0.22 (dotted), 0.30 (dash-dotted hne). 



The variable ^(a:,p|.), first suggested in [17], relates the (very different!) dependencies on longitudinal and transverse 
momenta so that the factorization of these dependencies can be only approximate. It is in some sense similar to the 
Nachtmann variable in DIS which controls some (kinematical) part of higher twist contributions. The reason for the 
apperance of such a variable is deeply related to the general properties of our model like Lorentz invariance [ITl [35] 
and the on-shellness of partons. 

Using as input for ff{x) the LO parameterization of [37] at the scale 4GeV^, we obtain for u and (i-quarks the 
results shown in Fig. [ij The lower part of this figure is shown again, in a different scale in Fig. [2| 
We make the following observations: 

i) For fixed x, the pT-distributions are similar to the Gauss Ansatz ff{x,pT) cx exp(— (pt))- This is an 
interesting result, since the Gaussian shape is supported by phenomenology [38 . 

ii) The width depends on x. This result reflects the fact, that in our approach the parameters x and pr are 
not independent due to rotational symmetry. 

in) The Figs, [ij [5] suggest that the typical values for transverse momenta, (pt) ^ O.OlGeV^ or {pt) ^ O.lGeV. 
These values correspond to the estimates based on the different analyses of the structure function F2(x,(3^) [24j. 
On the other hand, much larger values ~ OAGeV^ are inferred from SIDIS data referring to comparable scales 
[38'-'40^. Note also that in the statistical model of TMDs [36*, '4r the parameter (pt) may be interpreted as an effective 
temperature of the partonic "ensemble" It is instructive to compare this number to the lattice calculations [43] 
of the QCD phase transition temperature T ^ 175 MeV. 
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III. THE POLARIZED TMDS 



All polarized leading-twist T-even TMDs are described in terms of the same polarized covariant 3D distribution 
H(jP). This follows from the compliance of the approach with relations following from QCD equations of motion [23] , 
As a consequence all polarized TMDs can be expressed in terms a single "generating function" K^{x^y^t) as follows 



^i(^,Pt) 
^?(^,Pt) 



2x (( 



1 / 

— X 

X \ 

-i' 

X \ 
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m \ 

m) 

m 
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X K^{x,pt) 
X K^x.pt) 
X K^x.pt) 
X K^x.pt) 
X K^x.pt) 



(11) 



with the "generating function" K^{x^pt) defined (in the compact notation introduced in Eq. (32) of [23^) by 



K1{x,pt) = M^x J d{/} , d{/} 



dpi g^p") j fp' 

p° p° + m \ M 



P 



(12) 



We remark that in order to rewrite g\{x^pT) = / d{p^} [{p^ + m)xM — p^] [23^ as shown in the first equation of 
(11), we used the identity = (p|. + x^M^ + m?) / {2xM) valid under the p^-integral, which holds because in the 
model the partons are on-shell, i.e. Pq— p\ — = m^, and p^ — p^ = xM due to the delta- function. The expressions 
for the other TMDs in (11) can be read off directly from Eqs. (35-38) in [23 . Using the identity ([5| we perform the 
-integration in (12) and obtain for the generating function, which depends on x, p^ only via p^ defined in Eq. 



K\x,pt)=M^- 



m{p^) 



m 



f 



- xM 
2 



(13) 



From (11 13) it is clear that we can predict all polarized TMDs if we know H^{p^). The polarized 3D momentum 
distribution H^{p^) could be determined in principle from any polarized TMD, but the helicity parton distribution 
function g^x) plays a special role, because its x-dependence is known. The connection of gl{x) and H^{p^) was 
derived previously in [18l|2T]. In order to make this work self-contained we present here an independent derivation. 

We start from the expression for gl{x) which follows from (11 ) and proceed as in Eq. (|6|, i.e. we neglect m and use 
spherical coordinates such that p^ = p cos and p|. = p'^ sin ^, i.e. 

d^P ^Ta/^^/^2,^2 ^2 „:„2 m r /^^ " ^ COS l9 

2p2 



M 

{x^M^-p^ sin' 0)6 



(14) 



= ^ d0 y H\p) J dcOsO 

= 27rM'^ j dpH\p)^{xM -p)g(^p-^xM 



p — p cos 
M 



where the 6-function emerges in the same way it did in Eq. (p|. Consequently, we obtain 



gl{x) = 27rM'x 



M 



xM 



V 



dp 
p 



m{p) 



dpm{p) 

\xM 

Differentiating this equation and integrating it by parts gives the following results 



(15) 



dq^(x) ^ 
X ^ = 271 M^x 



dx 



gl{y) = 2^M^x 



2xM j ^H\p)- J dpH'^{p)-^H' 

\ 
J 




M 
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dp 
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dpH\p) 



\xM 



\xM 



(16) 
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FIG. 3: The TMD 5'?(x,Pt) for u- {upper panel) and d-quarks {lower panel). Left panel: Pi(x,pt) as function of x for 
Pt/M — 0.10 (dashed), 0.13 (dotted), 0.20 (dash-dotted hne). The sohd hne corresponds to the input distribution g\{x). 
Right panel: ^i(x,Pt) as function oipr/M for x = 0.15 (soUd), 0.18 (dashed), 0.22 (dotted), 0.30 (dash-dotted hne). 



Now we see that if we take the linear combination 2 dy gliy) / y -\- 3 gi{x) — x gWx) the integral terms from the last 



three equations cancel out, and we obtain for p^H^{p) at p = the expression 



^x^M^mi^x 



y 9i{y) + 3^J(x) - X 



which confirms previous works [T8l[2T]. For the generating function ( |13| ) we obtain, in the limit m 

2 



K^{x,^t) 



2 ? 



« y 



9l{y) + i9m-^ 



1 



Pr 



(17) 
0, the result 

(18) 



and from ( 11 ) we obtain (in agreement with the result reported in the proceeding [25 which was derived independently) 



91{x,tpt) 



(19) 



In (18 



and also below in (|20|) we use the variable ^ = ^(x,pt) as defined in (10). 



Eq. (19) yields for ^i(x,pt), with the LO parameterization of [44 for g^x) at 4GeV^, the results shown in Fig. [s] 



The remarkable observation is that gl{x^pT) changes sign at the point pt = Mx, which is due to the prefactor 

2^ 



2x-^ = x{l 



f Pt Y 
\MxJ 



-2p^/M 



(20) 



in (19). The expression in (20) is proportional to the quark longitudinal momentum p^ in the proton rest frame, 
which is determined by x and pt, see Eq. This means, that the sign of gi{x^pT) is controlled by sign of p^. 

To observe these dramatic sign changes one may look for multi-hadron jet-like final states in SIDIS. Performing the 
cutoff for transverse momenta from below and from above, respectively, should effect the sign of asymmetry. 
There is some similarity to g2{x) which also changes sign, and is given in the model by the expression [21] 



m 



-P 
M 



po ■ 



(21) 



The (5— function implies that, for our choice of the light-cone direction, large x are correlated with large and negative 
j9^, while low X are correlated with large and positive p^. Thus, g2{x) changes sign, because the integrand in (21) 
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FIG. 4: The TMDs hl{x, pt), (/^^^(x, px), hf:}{x,pT) for u- and d-quarks. Left panel: The TMDs as functions of x for 
Pt/M = 0.10 (dashed), 0.13 (dotted), 0.20(dash-dotted Unes). Right panel: The TMDs as functions oipr/M for x = 0.15 
(sohd), 0.18 (dashed), 0.22 (dotted), 0.30 (dash-dotted hues). 
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changes sign between the extreme values of p^. Let us remark, that the calculation of g2{x) based on the relation (21 ) 
well agrees [19 with the experimental data. 

The other TMDs (11) can be calculated similarly and differ, in the limit m ^ 0, by simple x-dependent prefactors 



K\x,pt)^ 
-K^{x,pt)^ 



(22) 



1 



K\x,pt) 



The resulting plots are shown in Fig. Ul We do not plot hf^ since this TMD is equal to — in our approach [23] . 
Let us remark, that g^x^ pr) is the only TMD which can change sign. The other TMDs have all definite signs, which 
follows from hi, 22). Note also that pretzelosity /i^^(x,pt), due to the prefactor 1/x, has the largest absolute value 
among all TMDs. 



IV. CONCLUDING REMARKS 

We have studied relations between the TMDs (x, Pt), Pt), hl{x, Pt), Pt), P^)' ^it(^' P^) 

the PDFs fl{x)^g\{x). These relations follow in the covariant parton model as a consequence of Lorentz invariance 
and rotational symmetry of parton momenta in the nucleon rest frame. These relations used with parameterizations 
of PDFs give numerical predictions for all T-even (in this work: leading twist) TMDs. In fact we have demonstrated, 
that in the framework of the model there are fundamental distributions covariant momentum distributions G^{p) and 
H^{p) from which all the considered unpolarized and polarized PDFs and TMDs can be uniquely obtained. And vice 
versa, from any PDF or TMD one can unambiguously calculate the corresponding or H^. 

Some of our results are compatible with the results of the recent paper [35 . In spite of some differences, both 
approaches have an important common basis consisting in the Lorentz invariance. For a more detailed comparison of 
the two approaches we refer to [24] . Our predictions are consistent also with the results obtained in the recent study 
[36] . some quantitative differences between these two approaches are discussed in the cited paper. 

To conclude, let us remark that an experimental check of the predicted TMDs requires care. In fact, TMDs are 
not directly measurable quantities unlike structure functions. What one can measure for instance in semi-inclusive 
DIS is a convolution with a quark fragmentation function. This naturally "dilutes" the effects of TMDs, and makes it 
difficult to observe for instance the prominent sign change in the helicity distribution, see Fig. [3| A dedicated study 
of the phenomenological implications of our results is in progress. 
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